Optimal Real-time Spectrum Sharing between 
Cooperative Relay and Ad-hoc Networks 

Yin Sunt, Student Member, IEEE, Xiaofeng Zhongt, Member, IEEE, Tsung-Hui Chang*, Member, IEEE, Shidong 
Zhou*, Member, IEEE, Jing Wang*, Member, IEEE, and Chong-Yung Chi*, Senior Member, IEEE. 



o 

<N 

a 



> 

i— H 

CO 
00 

o 

O 



Abstract — Optimization based spectrum sharing strategies 
have been widely studied. However, these strategies usually 
require a great amount of real-time computation and significant 
signaling delay, and thus are hard to be fulfilled in practical 
scenarios. This paper investigates optimal real-time spectrum 
sharing between a cooperative relay network (CRN) and a nearby 
ad-hoc network. Specifically, we optimize the spectrum access 
and resource allocation strategies of the CRN so that the average 
traffic collision time between the two networks can be minimized 
while maintaining a required throughput for the CRN. The devel- 
opment is first for a frame-level setting, and then is extended to 
an ergodic setting. For the latter setting, we propose an appealing 
optimal real-time spectrum sharing strategy via Lagrangian dual 
optimization. The proposed method only involves a small amount 
of real-time computation and negligible control delay, and thus 
is suitable for practical implementations. Simulation results are 
presented to demonstrate the efficiency of the proposed strategies. 

Index Terms — Cognitive radio, relay network, ad-hoc network, 
resource allocation, real-time control, spectrum sharing, collision 
prediction. 



I. Introduction 

In recent years, spectrum sharing between heterogeneous 
wireless networks has been recognized as a crucial technol- 
ogy for improving spectrum efficiency [1], [2] and network 
capacity [3], [4]. There are two major models for spectrum 
sharing presently, namely, the open sharing model and the 
hierarchical access model [1], [2]. In the open sharing model, 
each network has equal right to access the same spectrum 
band, e.g., the unlicensed band, and there is no strict constraint 
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on the interference level from one network to its neighbors. 
In the hierarchical access model which consists of a primary 
network and a secondary network, the secondary network, i.e., 
cognitive radio, dynamically accesses the spectrum provided 
that the primary users' transmission is almost not affected [1]. 
In either model, the inter-network interference make spectrum 
sharing a challenging task, especially when there is no explicit 
coordination between the coexisting networks. 

To address this interference issue, cognitive spectrum access 
strategies have been proposed [5]— [12] for the hierarchical 
access model. While these works focus on MAC-layer spec- 
trum access, there have been works focusing on physical- 
layer resource allocation of secondary networks, where strict 
constraints are imposed to limit the induced interference to the 
primary users; see [13]— [15] and also [16] for an overview. 
Joint optimization of spectrum access and resource allocation 
was studied in [17] for an open sharing model that considers 
spectrum sharing between an uplink orthogonal frequency- 
division multiplexing (OFDM) system and an ad-hoc network. 

However, some implementation issues were rarely consid- 
ered in these optimization based spectrum sharing strategies. 
First, real-time optimization can be computationally quite 
demanding for realistic wireless networks [18]. Second, spec- 
trum sensing and channel estimation are usually performed at 
spatially separate nodes, which requires information exchange 
between these nodes before solving the optimization problem. 
The resultant computation and signaling procedure usually 
lead to significant control delay, making these strategies hardly 
be fulfilled in practical scenarios. Therefore, spectrum sharing 
strategies with little real-time calculation and small signaling 
delay are of great importance for practical applications. 

In this paper, we study spectrum sharing between an uplink 
broadband cooperative relay network (CRN) and an ad-hoc 
network, as illustrated in Fig. 1. The CRN is composed of 
a source node (e.g., a mobile terminal), a relay node, and a 
distant destination node (e.g., a base station (BS)). The CRN 
adopts a two-phase transmission protocol for each frame: in 
the first phase, the source broadcasts an information message 
to the relay and BS; in the second phase, the relay employs 
a broadband decode-and-forward (DF) strategy [19]-[22] to 
forward the message. The source will transmit a new message 
to the BS in the second phase as well. In order to communicate 
with the distant BS, the source and relay transmit signals with 
peak powers, which, however, will induce strong interference 
to nearby ad-hoc links operating over the same spectrum 
band. The ad-hoc transmitters, e.g., wireless sensor nodes, 
have relatively low transmission powers due to their short 
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Fig. 1. System model. In Phase 1 (left plot), the source node broadcasts information to the relay and destination, and the transmitted signal 
interferes with the nearby ad-hoc links; in Phase 2 (right plot), the source and relay transmit signals to the destination simultaneously and 
both transmitted signals interfere with the ad-hoc links. 



communication ranges, and thus their interference to the relay 
and destination can be treated as noise. Such an asymmetrical 
interference scenario is known as the "near-far effect" [3] . Our 
contributions in this paper are summarized as follows: 

1) We first consider the case that the source and relay nodes 
perform spectrum sensing at the start of each frame. By 
modeling the ad-hoc traffic in each ad hoc band as inde- 
pendent binary continuous-time Markov chains (CTMC) 
[23], the average traffic collision time between the two 
networks can be obtained based on collision prediction. 
We formulate a frame-level spectrum sharing design 
problem that jointly optimizes physical-layer resource 
allocation and MAC-layer cognitive spectrum access of 
the CRN such that the average traffic collision time is 
minimized, while maintaining a required throughput for 
the CRN. 

2) The formulated spectrum sharing problem is a difficult 
nonconvex optimization problem with no closed-form 
expression for the objective function. We first derive 
the optimal spectrum access strategy, based on which, 
we show that the resource allocation problem can be 
reformulated as a convex optimization problem. To solve 
this problem in a low-complexity manner, we present a 
Lagrangian dual optimization method, which has a linear 
complexity with respect to the number of sub-channels. 

3) As common issues of existing frame-level transmission 
control strategies, the developed frame-level spectrum 
sharing strategy requires excessive computation and sig- 
naling process, which may cause considerable control de- 
lay and is not suitable for real-time implementations. To 
overcome these issues, we further formulate an ergodic 
spectrum sharing design problem based on a long term 
average CRN achievable rate and a long term average 
traffic collision time. By exploiting the Lagrangian dual 
optimization solution, we develop a low-complexity, real- 
time implementation method for obtaining the optimal 
spectrum sharing strategy. The proposed strategy is ap- 
pealing because most computation tasks are accomplished 
off-line, leaving only simple tasks for real-time com- 
putations. In addition, the computation and signaling 



procedure are carefully designed to maximally reduce 
the control delay. This method can accommodate an 
additional spectrum sensing in Phase 2 of each frame 
to further improve the accuracy of collision prediction. 

The spectrum sharing strategies proposed in this paper 
differs from that reported in [17] in three aspects: First, our 
interference metric is more practical in the considered strong 
interference scenarios (see Remark 1 in Section III-A for more 
details). Second, the spectrum sharing design problem of our 
DF based CRN is more difficult compared to that of point- 
to-point uplink system considered in [17]. Finally, an optimal 
real-time implementation method is proposed for the ergodic 
spectrum sharing design strategy, which is never reported in 
the literature before. 

The proposed strategies may provide potential spectrum 
sharing solutions for various application scenarios; e.g., the 
coexistence between unlicensed WiMAX (IEEE 802.16) and 
WiFi (IEEE 802.11) [24], the coexistence between a relay- 
assisted cellular network and ad-hoc networks including mo- 
bile ad-hoc networks [18] and peer-to-peer communication 
networks [25], and, in the domain of military communications, 
the coexistence between the broadband tactical backbone 
network and local ad-hoc networks, such as sensor networks 
and tactical mobile ad-hoc networks [26]. 

For ease of later use, let us define the following notations: 
The probability of event A is denoted by Pr{A}, and the 
probability of event A conditioned on event B is denoted 
by Pr{A|£?}. The indicator function of event A is given 
by 1(A). E^jX} represents expectation of X over random 
variable u, and E{X|Z = z} is the conditional expectation 
of X given Z = z. We denote || x ||2 as the Euclidean 
norm of vector cc, and denote tt(S) as the size (measure) of 
a set, e.g., 7r([a, b]) = b — a. The projections of x on the 
sets [0, oo ) and [0,y] are denoted by = max{x,0} and 
[x]q = min{max{x, 0}, y}, respectively. 

The remaining parts of this paper are organized as follows: 
In Section II, the system description and the formulation of 
frame-level spectrum sharing problem are presented. Section 
III presents a Lagrangian dual optimization method to resolve 
the frame-level spectrum sharing problem. The ergodic spec- 
trum sharing problem and its real-time implementation method 



are presented in Section IV. Conclusions are drawn in Section node transmits over the nth sub-channel in Phase 1, and 
V. 1^ C [aTf,Tf] denote the time set that the source and 

relay nodes transmit over the nth sub-channel in Phase 2, 
II. System Description and Problem Formulation for n = 1, . . . , N. As the example in Fig. 2 shows, and 



A. Signal Model 

We assume that the CRN employs a broadband multi-carrier 
air interface, where all nodes transmit and receive signals over 
TV parallel sub-channels, denoted by M = {1,2,..., TV}. The 
term "sub-channel" here represents either a frequency subband 
or a group of consecutive subcarriers in OFDM systems [1]. 
The ad-hoc links operate in M non-overlapped frequency 
bands, denoted by M = {1,...,M}. Moreover, the mth 
ad-hoc band overlaps with the sub-channels of the CRN in 
the set Mm ^ M for m = 1, . . . , M, where M m satisfies 
Um=i Nm = M and M p f]M q = for p ^ q. An example 
with M = 2 and N = 6 is illustrated in Fig. 2. 

In the temporal domain, all nodes in the CRN transmit 
and receive signals in a frame-by-frame manner, where each 
frame has a fixed duration Tf. The source node performs 
spectrum sensing at the start of each frame to detect the 
ACTIVE/IDLE state of each ad-hoc band. Perfect sensing and 
negligible sensing duration is assumed in this paper, which is 
reasonable for moderate sensing signal-to-noise ratio (SNR) 
and common value of the frame duration Tf l . Since the 
computation capability of the source node is quite limited, the 
source node cannot determine the optimal spectrum sharing 
design. Therefore, it forwards the obtained sensing results to 
the destination. The destination node computes the optimal 
spectrum access and resource allocation strategies of the CRN, 
and feeds them back to the source and relay nodes. As one may 
have noticed, the computation and signaling procedure may 
cause considerable control delay between spectrum sensing 
and data transmission. In this work, we will explicitly take 
into account this control delay in the spectrum sharing design. 
As illustrated in Fig. 2, we use STf to denote this control 
delay, where S G [0, 1). For large value of S, the transmission 
time of the CRN (1 — S)Tf becomes quite small, leading to a 
significant performance degradation. In Section IV, a real-time 
spectrum sharing strategy with negligible control delay will be 
presented. 

In practice, the relay node operates in a half-duplex mode 
[19]. Therefore, each frame consists of 2 phases: In Phase 1, 
the source transmits signal to the relay and destination via 
a broadcast channel; in Phase 2, the source transmits a new 
information message, and, at the same time, the relay uses the 
DF relay strategy to forward the information message received 
in Phase 1 to the destination, which forms a multiple-access 
channel. These scenarios are illustrated in Fig. 1. The time 
durations of Phase 1 and Phase 2 are set to aTf and (1 — a)Tf, 
respectively, where a G (£, 1)- 

We assume that the source and relay nodes can switch 
on and off their transmissions freely over each sub-channel. 
Let iffl C [STf, aTf] denote the time set that the source 

lr rhis assumption was confirmed in [8] for a WLAN energy detector. When 
the sensing SNR is 5 dB, a sensing duration of less than 5/xs is sufficient to 
achieve a detection error probability of 1CT 5 . If we choose a frame duration 
from 500/^s to 2ms, the sensing duration is fairly small. 



(2) 

In each may be a union of several disjoint transmission 
time intervals. Note that the source node cannot transmit 
during [0, <5T/], owing to aforementioned control delay. For 
convenience, let us define 
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which represent the fractions of the CRN transmission time in 
Phase 1 and 2 of each frame, respectively. 

We assume that the wireless channels of source-relay, 
source-destination, and relay-destination links are block-fading 
[27], which means that the channel coefficients remain static 
within each frame, and can change from one frame to another. 
Let hfy denote the frequency response of sub-channel n 
between transmitter i and receiver j, where i G {s, r} and 
j G {r, d} (i ^ j), in which s, r, d stand for the source node, 
relay node and destination node, respectively. The interference 
plus noise at the relay and destination nodes are modeled 
as independent, zero mean, circularly symmetric complex 
Gaussian random variables, with N£ and N d denoting the 
respective peak power spectral densities (PSD) over the nth 
sub-channel (i.e., the weak interference from the ad-hoc net- 
work to the CRN is treated as noise). Hence, the quality of the 
wireless links can be characterized by the normalized channel 
power gains g n r 4 \h n r \ 2 /N n W, g n d 4 \h n d \ 2 /N n W, and 
9n d ~ \ h n d \ 2 / N n W > where w is the bandwidth of each 
sub-channel. For broadband DF CRN with TV parallel sub- 
channels, the following rate is achievable [19, Eq. (45)] 
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where Ps]l and Ps% denote the transmission powers of the 
source over sub-channel n in Phase 1 and Phase 2, respec- 
tively; and P r?n is the transmission power of the relay over 
sub-channel n in Phase 2. The achievable rate Rcrn in (2) is a 
concave function of the transmission power and time variables 
{P&,p£l,Pr,n,0£\ { n\neM} [28, p. 89]. 

The ad-hoc traffic over the mth band is modeled as an inde- 
pendent, stationary binary CTMC X m (t), where X m (t) = 1 
(X m (t) = 0) represents an ACTIVE (IDLE) state at time 
t. The holding periods of both ACTIVE and IDLE states 
are exponentially distributed with rate parameters A and /i, 
respectively. The probability transition matrix of the CTMC 
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Fig. 2. Time-frequency transmission structure of the CRN and the ad-hoc network. In this figure, the number of CRN sub-channels is 6 
(N = 6) and the number of ad-hoc bands is 2 (M = 2). 



model of band m is given by [23, p. 391] 



P(t) 



1 



A + /i 



H — jie 



-(A+/x)t 



A + /ie 
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where the element in the zth row and jth column of P(t) stands 
for the transition probability Pr{X m (t + r) = j — l|X m (r) = 
i — 1} for z,j G {1,2}. This CTMC model has been used 
in many theoretical spectrum sharing studies and verified by 
hardware tests; see [7]-[ll], [17]. In practice, the parameters 
A and /i can be estimated by monitoring the ad-hoc traffic in 
idle frames of the CRN [7]. 

B. Traffic Collision Prediction and Interference Metric 

We utilize the average traffic collision time between the 
CRN and the ad-hoc network as the metric of interference 
experienced by the ad-hoc links. Since the ad-hoc nodes are 
near the source and relay nodes, the ad-hoc links would suffer 
from communication errors, whenever the ad-hoc transmission 
happens to collide with the CRN traffic 2 . Let x m G {0, 1} 
denote the sensing outcome for the rath ad-hoc band, i.e., 
X m (0) = x m . Given x m , one can predict the average traffic 
collision time based on the CTMC model in (3). Specifically, 
the average traffic collision time over the rath ad-hoc band is 
given by 



l(X m (t) = l)dt 




l(X m (t) = l)dt 



X m (0)=x m },(4) 



where l(X m (i) = 1) is the indicator function of event 
X m (t) = 1. It is worthwhile to note that U n eAT m ^ in ( 4 ) 

2 In general, the transmission error probability of an ad-hoc band is an 
increasing function of the average traffic collision time. In particular, if the 
ad-hoc transmission is uncoded in the physical layer or the ad-hoc code 
block length is quite short, the transmission error probability caused by traffic 
collisions is approximately proportional to the average traffic collision time. 



is the time set that the CRN is transmitting in at least one 
sub-channel in Af m , where Mm is the set of sub-channels 
overlapping with the rath ad-hoc band. This reflects the fact 
that, in the considered strong interference scenario, the ad-hoc 
transmission in the rath band is disrupted, even if the CRN 
transmits in only one sub-channel of M m . 

Each of the expectation terms in (4) can be calculated as 
follows 
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X m (0) = x ri 
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J LJ t %) 
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/ Pr{X m (t) = l\X m (0)=x m }dt (5) 

J U ^ 



By (5), the total traffic collision time summed over all the 
ad-hoc bands is given by 



rn=l 



/ Pr{X m (t) = l|X m (0) = x m } dt 
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L nGM 



u 



Pv{X m (t) = l\X m (0) =x m }dt 

tt(2) 



•(6) 



C. Frame-level Spectrum Sharing Design 

The goal of the CRN is to optimize the source and relay's 
transmission powers Ps] n -,Ps 2 n and P r , n , and their spectrum 



access strategies and In J , such that the total traffic 
collision time in (6) is minimized, while a minimum uplink 
throughput i^min can be maintained. This joint spectrum ac- 
cess and resource allocation design problem can be formulated 
as the following optimization problem: 
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(Id) 



P£lP$%Pr,n>0, n=l,...,7V, (7e) 
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where Rcrn and / are given by (2) and (6), respectively, (7g) 
follows from (1), P^ax anc * ^m ax m (7c) and (7d) denote the 
power constraints at the source and relay nodes, respectively. 

III. Optimal Frame-level Spectrum Sharing 
Solution 

Problem (P) is difficult to solve because the objective 
function / has no closed-form expression. Fortunately, this 
issue can be resolved by analyzing the optimal transmission 
time sets in the two phases of each frame, i.e., I n 1 ^ and I n 2 \ 
from which Problem (P) can be reformulated as a convex 
optimization problem, as we will present in this section. A 
Lagrangian dual optimization method is also proposed to 
obtain an optimal solution of (P) efficiently. 

A. Reformulation of Problem (P) 

The key idea that makes this convex reformulation possible 
is to examine the optimal spectrum access and £ ] m 
Problem (P). In particular, it can be shown (in Lemma 1 
below) that the optimal spectrum access must satisfy the 
following two principles: 

1) In both Phase 1 and Phase 2, the source and relay nodes 
should transmit as soon (late) as possible if the sensing 
outcome is IDLE (ACTIVE); 

2) The CRN should have identical spectrum access strategy 
for the sub-channels overlapping with the same ad-hoc 
band; that is, 1^ = 1^ for all p : q G jV m , where m G M 
and i G {1,2}. 

Principle 2) shares the same flavor of interference alignment 
technique in [29] since both of them align the transmissions 
to reduce the interference to the ad-hoc links. Let us define 



W A 



max j 0^,71 G Mnj 



(8) 



as the largest transmission time fraction over the sub-channels 
Mm in phase i. The optimal spectrum access strategy is given 
by the following lemma, which is proved in Appendix A. 



Lemma 1 For any given transmission time fractions {9$ G 
[0,a — 6], 9$ G [0, 1 — a}}^ =1 and transmission power 
{Pi,n, Ps 2 l-> P r , n }n=i * n Problem (P), we have that: 

1 ) The optimal spectrum access strategy in Phase 1 is given 
by = [5T f , (5 + §$)T f ] for all n G M m if the 
sensing outcome is x m = 0, and is given by iffl = [(cm — 
9$)Tf,aTf] for all n G Af m if the sensing outcome is 

2 ) The optimal spectrum access strategy in Phase 2 is given 
by ll 2) = [aTf, (a + 9$)T f ] for all n G M m if the 
sensing outcome is x m = 0, and is given by In = [(1 — 
9^)Tf^ Tf] for all n G Mm if the sensing outcome is 

An example that illustrates the spectrum access strategy of 
Lemma 1 is shown in Fig. 3. 

Remark 1: A similar result of Lemma 1 was re- 
ported in [17] for a different interference metric, which 
cumulates the traffic collisions for all the sub-channels 
overlapping with the same ad-hoc band. For example, 
EnejV m E { x m(t)dt X m (0) = x m | represents the av- 
erage traffic collision time over the rath ad-hoc band in 
[17] instead of our interference metric given in (5). Our 
interference metric is more practical than the one reported 
in [17], because, in strong interference scenario, the ad-hoc 
transmission in the rath band is disrupted, no matter that the 
CRN is transmitting in either one or more sub-channels of 
Mm- Since transmitting in more sub-channels of M m will not 
further worsen the interference, the CRN should transmit over 
all the sub-channels of M m simultaneously to increase data 
rate, which is different from the result of [17]. 

According to Lemma 1 , the integration region in each term 
of (6) is a simple time interval. In order to simplify (6), we 
define for G [0,a-S\, 

(j)(l)(9\Xm 0) 

= / Pr(X m (t) = l|X m (0) = 0)& 

J[6T fi (6+0)T f ] 
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Fig. 3. An example illustrating the spectrum access strategy of Lemma 1. The sensing outcomes of Frame 1 are x\ — 1 and X2 = 0, and 
the sensing outcomes of Frame 2 are x\ — and X2 = 1. 
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The interference metric in (6) can be simplified as 



T f . (12) 



M 

I= E [^(1) {^;x m )+^ (2) (e^;x m )] . (13) 



m=l 



It is easy to verify that the functions 4>n\(0;x) and 0( 2 ) (0; #) 
in (9)-(12) are strictly increasing and strictly convex functions 
of 0. Thus / in (13) is a convex function of 



(i) 



and 



0m . Further, the control delay STf degrades the interference 
mitigation performance of the spectrum access strategy in 
Lemma 1, as 0m(0;O) is strictly increasing in S. 

On the other hand, it follows from (2), (8) and Lemma 1 
that the constraint (7b) of (P) can be equivalently expressed 
as 



By (13) and (14), Problem (P) is equivalent to the following 
problem: 
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which is a convex optimization problem. While problem (17) 
can be solved by interior-point methods, we present in the next 
subsection a low-complexity Lagrangian dual optimization 
(14) method. 
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+fl®log 2 1- 



p(2) ~s,d 
r s,ny n 



p(l) ^s,d 
^Pr,ng n d 



5(2) 



Suppose that problem (17) is strictly feasible. Then, accord- 
(15) ing to the Slater's condition [28], the strong duality holds for 
(17). Hence we can alternatively consider the following dual 
optimization problem 



(16) 



max 



mm 



(18) 
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where V 4 {(P^ , p$ , P r , n , 0« , §$ ) 



5(2) 



L = 



>(i) 



>(2) 



0<^ } <a- 



^(i) 



< 1-a, P^,Ps,n,P r ,n > 0,nGjV,mGM} 

M 

+" 0(2) (^m^m) 



M 
m=l 



■5, < We first solve (20) to obtain the optimal ratio -0y . Specifi- 
cally, if Pin > 0, then equality in (20) holds, and an optimal 

p(i) 

is equal to the positive root x of the following quadratic 



w 



a 



(-Rmin — Rl ) + (-Rmin — R-2 ) 



N 



/ j \ s,n^ ± s,nJ x rr 



+ 7? 



TV 



EP -P r 
J r,n 1 rr 



, (19) 



is the partial Lagrangian [30] of (17), (> 0, a > 0, e > 0, and 
77 > are the dual variables associated with the constraints 
(17b), (17c), (17d), and (17e), respectively. As will be shown, 
the inner minimization problem of (18) has closed-form solu- 
tions for Pin, Pi n and P r n for n = 1, . . . , N, and 0$ and 
6m for m = 1, . . . , M. Hence, the computational complexity 
for solving the inner problem is only linear with respect to N 
and M. Moreover, the outer maximization problem of (18) 
only involves four optimization variables (£, cr, £,77), which 
is much smaller than the number of variables of the primal 
problem (17). 

Suppose that a dual variable v = (£, a, e, rf) T is given. Let 
us present the closed-form solutions of the inner minimization 
problem of (18). Because the inner problem is convex, the 

optimal (Pi,n, Psfn, Pr,n,dm\ Qm*) for fixed dual variable v 
must satisfy the Karush-Kuhn-Tucker (KKT) conditions [28] 
of the inner problem, which can be expressed as: 



dL 

U± s n 



dL 

dp^l 



dL 

dPr.n 



dL 



Cmax{^' r ,^ d } 



l+max{^^n' d }fif )ln2 (l+g 



,dPl 



In 2 



(g s n d 



J = 0ifP$>0 , om 

i>oif^=o (20) 

&9n 



+ E 



= if Pi 2 2>0 
> if Pi 2 2=0 



(21) 



V9n 



-1 1 s. 
1 +#n 



■ „r,d P r ,n 
-T yn ~ (2 ) 



ln2 



<9<9 



5(1) 



/( maX {^ 



P 

9n I 



(i) 
0(1) 



+ 



#0(1) JXmj 



< if 



= if Pr,n > 
> if Pr,n = 

p(l) 

J\ 9n 



(22) 



a(i) 



if ^ } G(0,a-(5) (23) 



> if 



2(i) 







dL 

»fl( 2 ) 



p(2)\ / p(2) 

l-Ls,n \ X ^ /»/ s.d-ts. 



> \ ^ /»/ s.d-fs^n \ \ ^ ,»/ s.d-fs^n r.d Pr,n 

neJV m \ ^ m / neJV™ \ ^ ^ 

< if = 1-a 
: if 0^ } e(O,l-a 



00(2) ^mj 



<5>6>. 



where f(x) = log 2 (1 + x) 



(l+x) In 2 



> if 6$ = 



in (23) and (24). 



equation 



ag, 



Cmax{^ r ,^ d } 
1 + max{^' r , gn d }x 1 + 0. 



: eln2. 



(25) 



If (20) has no positive root, then 



p(i) 



(i) 



0. Next, let us find the 



p(2) p 

optimal and 



-^y by solving (21) and (22). If P r?n > 

d(2) 
s , n 

^(27 

be obtained from (21) and (22) as 



the equality in (22) holds, and the optimal and can 



p(2) 

Urn 
P 

1 r,n 



C 



(e - V9n d /9n d ) In 2 5. 

p(2) 



1 



n + 



77 ln 2 



1 

r,d 



,r,d ^(2) ' 



where = max{x,0}. Instead, if P r?g 
from (21) and (22) that 



.(2) 



#2) 
P 

1 r,n 



(J 



£ln2 



1 



3(2) 



0, 



(26) 
(27) 

0, we obtain 

(28) 
(29) 



where (28) is actually the water-filling solution when the 
source directly communicates with the destination without the 
use of relay. 

The optimal 0^ and §ffl can be obtained by solving (23) 

p(2) 

s,n 

0^ 



and (24), respectively, provided that the optimal 



and 



Pr,n 
fl(2) 



have been obtained from (25)-(29). By substituting 

p(l) p(2) p 

^7, and into (23)-(24), and by the definitions of 
and 0(2) (^; x) in (9)-(12), we can obtain the optimal 
values of 6$ and (fm as follows: For x m = 0, we have 1 



0(D 




(A+M)7> 



vf[9 S n d - 



+C/(max{<- r ,^ d }^ 



(30) 



J o 



0(2) 



1 



(A+M)7> 



1- 



neJVrr 



' ( 



p(2) 
g(2) 



a 



(31) 



(24) where [x]q = min{max{x, 0}, 7/}, and for x m = 1, we have 

! For the notational simplicity in (30)-(33), we have extended the definition 
of the natural logarithm \n(x) to that with \n(x) = — oo for x G (— oo, 0] . 
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Algorithm 1 The proposed Lagrangian dual optimization 
algorithm for solving (P) 

1: Input system parameters (TV, M, a, P^ ax , P£ ax , T f , 
Rmin), the ad-hoc traffic parameters A, /i, the channel 
quality {^ ,r ,^' (i ,^ d }JL 1 , the sensing outcome 
{xm}%l =1 , the computation and signaling delay parameter 
S, and a solution accuracy e. 

2: Set the iteration number k = 1 ; initialize the dual variable 

3: Compute the optimal {pjj, Pi,n, ^r,n}^=i and 

{^ ) ,^ ) }m=i according to (25)-(33). 
4: Update the dual variable v k +i according to (34) and (35). 
5: If || i/fe || 2 < e, go to Step 6; otherwise, set k = fc+1 

and return to Step 3. 
6: Output the optimal primal solution {Pi,n , ^i,n , P r ,n}^=i 
and{£ ) ,^ ) }^f = i- The optimal spectrum access strategy 
=1 can be obtained by Lemma 1. 



{lM 2) }f 




By substituting (30)-(33) into (25)-(29), the optimal 
Ps]l,Ps 2 l,P r ,n can then be obtained. 

What remains for solving (18) is to optimize the dual 
variable v = (£, a, £, 7/) T for the outer maximization problem. 
In view of that the dual function of (17) (the optimal value 
of the inner problem of (18)) may not be differentiable 
[31], we consider to update v using the subgradient method 
[32]. Specifically, at the fcth iteration, the subgradient method 
updates v by [32] 



(34) 



where the subscript k denotes the iteration number, s k is the 
step size of the kth iteration, and h(v k ) is the subgradient of 
the dual function, which is given by [31] 



h(v k ) 



(Rn 
(Rn 



Ri)/W 

RD/w 



V AT / p (l)* 



where P 



(i)* 



Z_-m= 

p(2)* 



s,n J 



1 - 1 r,n 



(35) 



and P* n are the optimal solution of the 
inner minimization problem (18) at iteration k, and R\ and P£ 
are the corresponding rate values in (15) and (16), respectively. 
It has been shown that the subgradient updates in (34) converge 
to the optimal dual point i/* as k —> oo, provided that the step 
size Sk is chosen according to a diminishing step size rule 



TABLE I. Parameter setting for Fig. 4. 



Parameter 


N 


M 


AT; 




max 


-'max 


a 


S 


Value 


2 


2 


1 


1 


1 


1 


0.5 


0.1 


Parameter 


X! 


X 2 


s.d 


s.r 


r.d 
01 


02 


s.r 
02 


r.d 
02 


Value 





1 


0.4 


1.3 


1.3 


0.5 


1.4 


1.4 



[32]. The convergence speed of the subgradient method can be 
improved if one further considers the acceleration techniques 
in [31], [33], [34]. In Algorithm 1, we summarize the proposed 
Lagrangian dual optimization algorithm of (P). 

C Simulation Results 

In this subsection, we provide some simulation results to 
examine the performance of the proposed CRN spectrum 
sharing strategy in (P). The parameters used in the simulations 
are listed in Table I. Since N = M = 2, we simply set 
M = {l}andA^ 2 = {2}. 

We compare the proposed CRN spectrum sharing strategy 
with two degenerated strategies, namely, the relay -free strategy 
and the sensing-free strategy. Similar to the work in [17], the 
relay-free strategy only considers the direct uplink transmis- 
sion from the source to the destination. To implement this 
strategy, we simply set g^ r = g% d = for Problem (P). 

In the sensing-free strategy, the average traffic collision time 
in (6) becomes 



rn=l 



u 



Pv{X m (t) = l}dt 



Pi{X m (t) = l}dt 



due to the lack of sensing result X m (0) = x n 
reduces to 



(36) 



It further 



I 



M , 



m=l 



(37) 



by means of the interference alignment principle in Lemma 
1. Then, the sensing-free strategy is obtained by solving (17) 
with the objective function in (17a) replaced by (37). 

Figure 4(a) presents the performance comparison results 
of normalized average collision time (I /Tf) versus required 
uplink spectrum efficiency R m i n / (NW). We observe from this 
figure that for R m in/(NW) < 0.07 bits/s/Hz, the relay-free 
strategy exhibits comparable performance with the proposed 
strategy; whereas for R m in/(NW) > 0.07 bits/s/Hz, the 
proposed strategy yields a smaller average traffic collision 
time. The performance improvement is attributed to the DF 
relay techniques. The sensing-free strategy always generates 
more traffic collisions than the proposed strategy, because it 
does not utilize the spectrum sensing outcomes. 

Figure 4(b) displays the optimal transmission time frac- 
tions {OP^O™ }n=i of (P) versus required uplink spectrum 
efficiency R m in/(NW). The CRN only transmits over sub- 
channel 1 when R m in/(NW) < 0.42 bits/s/Hz, because the 
sensing outcomes are x\ = and £2 = 1. For R m in/ (NW) > 
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Fig. 4. Simulation results of the proposed frame-level transmission control strategy of the CRN with parameters given in Table I. (a) 
Normalized average collision time (7/7/) versus required uplink spectrum efficiency R m in/ (NW), (b) optimized transmission time fractions 
versus required uplink spectrum efficiency R m in/(NW). 



0.42 bits/s/Hz, the CRN starts to transmit over both sub- 
channels to achieve more stringent throughput requirement 
Rmin/(NW). The maximal value of 6^ is a - S = 0.4, 
because the computation and signaling delay is STf = 0.17/. 
On the other hand, it is interesting that 0^ never achieves its 
maximum 1 — a = 0.5, even at the largest feasible value of 

(2) 

Rmin/(NW). This is because the transmission time n ' only 
contributes to R 2 , but not to R lm When 0^ = 0^ = 0.4, 
R\ achieves its maximum and Rcrn = min{iJi, #2} is 
constrained by R\ , the CRN cannot increase Ri by increasing 

M 2 ) 



After obtaining v*, we need to compute the transmission 
parameters based on real-time spectrum sensing and chan- 
nel estimation results. These real-time computation tasks are 
fulfilled carefully in order to minimize the computation and 
signaling delay STf. 

Another benefit of this ergodic setting is that it allows one 
more spectrum sensing at the beginning of Phase 2 to improve 
the accuracy of collision prediction. This, however, cannot be 
exploited in the frame-level setting, because Problem (P) must 
be solved before this additional sensing in Phase 2 is carried 
out. 



IV. Ergodic Spectrum Sharing Design and 
Real-time Implementation 

In the previous section, the optimal spectrum sharing design 
is obtained in each frame. As we mentioned in Section I, such 
a frame-level spectrum sharing strategy may encounter several 
implementation issues. First, the CRN has to solve Problem 
(P) within every frame, which may be computationally too 
demanding for realistic wireless networks. Moreover, the desti- 
nation node has to collect the spectrum sensing outcome from 
the source and/or relay nodes in order to solve (P). After 
solving (P), the destination node has to send the solutions back 
to the source and relay nodes. The computation and signaling 
procedure may cause a considerable control delay STf, leaving 
very short time for data transmission. 

These issues intrigue us to investigate more practical spec- 
trum sharing strategies that allow real-time implementations. 
The key ideas are 1) to reduce the amount of real-time 
computations and 2) to decrease the computation and signaling 
delay. Our approach is to consider an ergodic resource alloca- 
tion problem. Recall from Section III-A that the frame-level 
design problem (17) can be solved by the Lagrangian dual 
optimization method. We will show next that for the ergodic 
spectrum sharing design problem, one can compute the optimal 
dual variable 1/* off-line, and only some simple tasks are left 
for real-time computation. 



A. Ergodic Spectrum Sharing Design Problem 

Suppose that both the source and relay nodes perform one 
extra spectrum sensing at the beginning of Phase 2. When 
there is no sensing error, the sensing outcome for the mth 
ad-hoc band is denoted as X m (aTf) = y m G {0,1} with 
m = 1,...,M. This additional sensing outcome can be 
utilized in our spectrum sharing design problem to reduce 
traffic collisions, at the cost of extra computation and signaling 
delay in Phase 2, as shown in Fig. 5. For notational simplicity, 
the duration of this extra computation and signaling delay is 
also assumed to be STf. 

Let us define a network state information (NSI) as 

<*> = {9n d i9 S n r i9n d , X rn,yrn, 71 G N 771 G M} , 

which includes both the channel estimation and spectrum 
sensing results over the two phases. In the ergodic setting, 
the NSI u: varies across frames. We assume that the channel 
fading gains and the ad-hoc traffic are stationary and ergodic; 
furthermore, their statistical distributions are known to the 
destination node prior to transmissions. 

Similar to (6), the average traffic collision time for a frame 
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Fig. 5. Time-frequency transmission structure by Lemma 2. Here the sensing outcomes of Frame 1 are xi 
and the the sensing outcomes of Frame 2 are x\ = 0, xi = 1, y± = 1 and y2 = 0. 



1, x 2 = 0, 2/1 = and 2/2 — 1, 



with NSI u> is determined by 



/(«) 



Pr{X m (t) = l|X m (0) = a; m }di 



m=l\ J 

/ Pr{X m (t) = l|X m (aT / )= 2/m }(it , (38) 

^ U / 



+ 



where In denotes the set of CRN transmission time 
over sub-channel n in phase i given the NSI u. Note from 
(38) that, in contrast to (6), the collision time in Phase 2 
now depends on the sensing outcome y m . Since the NSI is 
stationary and ergodic across the frames, the long term average 
traffic collision time can be obtained by taking the expectation 
of over the distribution of the NSI uj [23], i.e., 



til ' 



Pr{X m (t) = l\X m (0) = x m }dt 



+/ 



Pr{X m (t) = l|X m (aT / ) = 2/m }dt .(39) 



Let us define the transmission time fractions 9n\u>) as in 
(1), and define 

0« (w) 4 max { flW («) , n € A/™ } , (40) 

for m = 1, . . . , M and z = 1,2. It is not difficult to show 
that the optimal spectrum access strategies stated in Lemma 
1 also hold true for the case with two spectrum sensings in 
each frame: 

Lemma 2 For any given transmission time fractions 
{^(w) G [0,a-i],4 2) M e [0,l-a-<S]}*f =1 we have 
that: 

1) The optimal spectrum access strategy of Phase 1 is given 
by = [6T f ,(6 + 0&\u>))T f ] (l£\u>) = [(a- 

0m^ (uj))Tf, aTf]) for all n G Af m , if the sensing outcome 
of Phase 1 is x rn = (x m = 1); 



2) The optimal spectrum access strategy of Phase 2 is given 
byl^\uj) = [(a + 5)T f ,(a + 5+9%\u>))T f ] (1^ (w) = 
[(1 - §£\u>))Tf,T f ]) for all n G M m , if the sensing 
outcome of Phase 2 is ym = (y rn = 1). 



An example of Lemma 2 is illustrated in Fig. 5. 
Define the two functions 

4>(2)(0;y m = 0) 
= / PT(X m (t) = l\X m (aT f ) = 0)dt 

J[(a+5)T f ,(G+a+5)T f ] 



\ ( e -(Mp.)6T f 



-(Mp)eTf_ 1 



A+/i [ (A + 
4>(2)(0;y m = 1) 

/ Pr(X m (t) = l\X m (aT f ) = l)dt 

J[T f -GT f ,T f ] 



(41) 



A 



H e 



A+/i 1 A (\+fJL)T f 



e (\+n)9T f _ 1 



T f , (42) 



where 6 G [0, 1 — a]. In accordance with Lemma 2, / in (39) 
can be simplified as 



( M 

I = E U , Yl [^(1) (*m M; X ™) +^(2) (^M; frr.) 



, (43) 



where 0m (0, x m ) has been defined in (9) and (10) for x m = 
and Xm = 1, respectively. 

The achievable average rate of the multi-carrier CRN can 
be shown to be 



Rcrn = min{R 1 ,R 2 } , 



(44) 
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where 



P$(u),P rin (u>) > 0,n e AT,meMy and 



(45) 



R 2 = WJ2 J2 E - 



id) 



(w)log 2 1 + 



+C j Miog 2 i + 



(46) 



and Ps]l{Lo),Ps%(w),P r , n {uj) are the transmission powers 
for a given NSI u. We should point out that the average rate 
Rcrn is not the ergodic data rate in the Shannon sense, 
but one achieved by taking the average over many adaptive 
channel coding blocks. The relay node will not transmit but 
queue up its received data from the source node [19], [35], 
until the channel quality and sensing outcome in Phase 2 is 
favorable. 

It follows from (43) and (44) that the ergodic spectrum 
access and resource allocation problem is 



{M 
ki)fe } M; *m) +0 ( 2)& 2) (<*>); \ 
m^l 

P r , n (u,),0d)(u,), 



(47a) 
(47b) 

(47c) 
(47d) 



S.t. R\ > P m in, ^2 > ^min 

U=i 

P^H,Pi 2 n )H,P r , n M >0, n = l,...,7V 

(47e) 

< < a-6,0 < 0${u) < l-a-5, 

m = l,...,M. (47f) 



5. Solving Problem (47) w Real-time 

In the sequel, we provide a real-time method to solve 
Problem (47) following the idea in Section III-B. Specifically, 
we solve the following dual optimization problem 



max min _L (48) 

C,<7,£,r7>0 V (pj 1 ? )(u,),P s ( 2 ? )(u;),P r;ri (u;),^)( U7 )^(2)( U7 )) G v / 



where V 4 { (pj^u;), Pi, 2 i(u;), P r ,»(w), ^(o;),^^) 
< ^(w) < a - 5, < ^ } (u>) < 1 - a - 5,P&(u) : 



1 { M 



+ ^(Pmin-Pl) + p^(Pmin #2) 



n=l 

N 



(49) 



is the partial Lagrangian of (47). 

Let us define u> £ 4 {^), <£<^), Xm (£), y m (£), 

n G A/", m G M} as the NSI of the M frame. Problem (48) is 
in general non-causal, because it requires the NSI realizations 
of future frames to solve the inner minimization 
problem. However, by making use of the statistical distribution 
of the NSI, we can solve Problem (48) in two steps: First, we 
optimize the dual variables v = (£, a, e, rf) T off-line based on 
only the statistical distribution of the NSI. Then, the primal 
solution is updated on-line according to the current NSI uj£. 
Most computation tasks are accomplished in the off-line dual 
optimization step, leaving only simple computations for real- 
time primal solution update, as detailed in the subsequent two 
subsections. 

1) Off-line Dual Optimization: Given a dual variable v = 
(C, 77) T , the optimal solution to the inner minimization 
problem of (48) for the NSI realization uj is provided as 
follows: 



The optimal 



and 



3^7 



of the inner 



minimization problem of (48) can be exactly obtained by (25) 



(29), with P s ( %,Pr,n and 6% replaced by P^(w),P r>n (a;) 
and §m respectively. The optimal 0$ (a;) can be obtained 
by either (30) or (32), depending on the sensing result x m . 
By (41) and (42), the optimal 0% \u>) is given as follows: If 

Vm = 0, 



1 



(\+fi)Tf 



3(2) 



otherwise, for 



= 1. 




J 



^ 2) (^) = 



{\+H)T f 



In- 



riEAfrr, 



9n 



I Pr,n(u) 

My 



5(2) 




By substituting (30), (32), (50) and (51) into (25)-(29), the 
optimal values of Ps}2(w), Ps 2 n(w) and P r>n (u;) can then be 
obtained. 

The optimal dual variable v* is obtained by a series of the 
subgradient update in (34), where the subgradient h(yk) is 



determined by 



h{v k )-- 



(Rn 

(Rn 

N 



Ri)/w 
R* 2 )/W 



N 



E E u feH}-F n 

n=l 



, (52) 



where Pi,n *(<*>), ^i,n and P* n (<jj) are the optimal so- 
lution of the inner minimization problem (48) for given dual 
variable and NSI realization uj, and R 1 and i? 2 are the 
corresponding rate values in (45) and (46), respectively. 

In the off-line dual optimization procedure, the true NSI 
realizations {u^}^ are not available. However, the statistical 
distribution of the NSI u is available. We can still compute 
the expectation terms in (45), (46) and (52) by means of 
Monte Carlo simulations. In particular, one may randomly 
generate a set of NSI realizations u> following the distribution 
of the channel fading and sensing outcome. Then, the inner 
minimization problem of (48) is solved for each artificially 
generated NSI realization ui. The expectation terms in (45), 
(46) and (52) can be obtained by averaging over these NSI 
realizations. Then, we can use the subgradient method in (34) 
to derive the optimal dual variable v* off-line. 

2) On-line Primal Solution Update: After v* is obtained, 
one can compute the optimal transmission parameters of Prob- 
lem (47) based on the true NSI {u^}^. More specifically, 
in Frame £, the CRN needs to solve the inner minimization 
problem of (48) for the optimal dual variable v = i/* and 
real-time NSI uj = uj £ . The computational complexity of this 
step is much smaller than that of off-line dual optimization, 
because each subgradient update involves solving the inner 
minimization problem of (48) many times for each of the 
artificially generated NSI realizations. On the other hand, spec- 
trum sensing and channel estimation are usually performed 
at spatially separated nodes, and thus additional information 
exchanges are required to accomplish the computation task. 

In the sequel, we present how to update real-time primal 
solution such that the computation and signaling delay STf 
can be substantially decreased. In practice, the BS is able to 
acquire the channel gain te' r (f),^' d (f)j^' d (<)}JLi through 
channel prediction before Frame £ starts, if the wireless 
channel varies slowly across the frames [36]. Given the 
channel gain {gi^ r (£) , gi^ d (£) , g% d (£)} , the BS can compute 



the ratios °\V , *<Z K and 



according to (25)- 



(29) in Frame i — 1, and then feed the results back to the 
source and relay nodes before the start time of Frame £. 
Once the source and relay nodes obtain the sensing outcomes 
x m (£) and y m (£), they can compute 0m\uie) and Om\ue) 
according to the closed-form solutions (30), (32), (50), and 

(51). By providing the source and relay nodes with - 



and p ^(^) ahead of time, the signaling procedure 



0m 5 " Ina e^\ wt 
will not cause additional control delay. As 0$(ujA and 
Om (ue) has closed-form solutions, the time delay STf for 
computing 0$ (cc^) and 6$ ] (u>e) is quite short. Therefore, the 



total computation and signaling delay is quite short. In each 
frame, the destination node needs to send 3A^ parameters to the 
source and relay nodes. Note that this amount of information 
exchanges are equal to that of conventional OFDM CRN 
without spectrum sharing [19], where the BS also needs to 

feed the ratios ^J^, and %^ back to the 

source and relay nodes. 

In practice, the source and relay nodes may not be able 
to compute 0$ (u £ ) and 0$ (u £ ) in (30), (32), (50), and 
(51), owing to hardware limitations. In this case, the BS can 
compute two possible values of §m\<*>e) in (30) and (32) in 
advance for both of the sensing outcomes x m (£) = 1 and 
x m (f) = 0. Similarly, m (^i) in (50) and (51) can also 
be computed in advance at the BS for both y m (£) = and 

y m (£) = 1. The BS sends #^ ™* 



and the two possible values of 0m\ui) and 9m\wi) to the 
source and relay nodes. After spectrum sensing, the source and 
relay nodes simply select the values of §m\ue) and §m\ue) 
according to the sensing outcomes x m (£) and y m (£), respec- 
tively. This strategy has two benefits: First, all computations 
are carried out at the BS, and source and relay nodes only 
need to select appropriate transmission parameters. Second, 
the source and relay nodes can start transmission right after 
spectrum sensing, which means the computation and signaling 
delay STf is negligible. In each frame, the destination needs 
to send 3Af + 2M parameters to the source and relay nodes, 
which is slightly larger than the strategy introduced in last 
paragraph. 

Remark 2: With the optimal dual solution v* obtained ahead 
of time, the transmission parameters of the CRN in Phase 
1 is determined solely by x m (£) but not y m (£). Hence, the 
primal solution to the ergodic transmission control problem 
can be computed in a causal manner. However, y m (£) cannot 
be exploited in the frame-level transmission control problem 
(P), because neither the dual optimal solution v* nor the 
primal optimal solution to (P) can be obtained in advance 
without y m (£). 

Remark 3: If the source transmits to the relay node in Phase 
1, and the sensing outcome becomes unfavorable in Phase 2, 
the relay node will queue up its received data from the source 
node, and wait for better transmission opportunity [19], [35]. 
In other words, data queuing at the relay node allows the CRN 
to exploit transmission opportunities in the two phases of each 
frame separately. 

Remark 4: In practice, sensing error may occur at the source 
and relay nodes, which means that the ACTIVE (IDLE) ad- 
hoc traffic state is mistakenly detected as IDLE (ACTIVE). 
In this case, the source and relay nodes may transmit in 
different time intervals, leading to extra collisions to the ad- 
hoc traffics. In spite of no transmission synchronization in 
source and relay nodes, the destination node is still able to 
decode the messages sent from the source and relay nodes. The 
destination node first decodes the message from the relay node, 
and then decodes the source's message, by means of sequential 
interference cancelation (SIC) decoding [19]. The performance 
degradation of the interference metric (i.e., long term average 
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(c) R m in/(NW) = 1.7 bits/s/Hz 



(d) R m in/(NW) = 2.8 bits/s/Hz 



Fig. 6. Simulation results of the proposed ergodic transmission control strategy of the CRN. (a) Normalized long term average collision time 
(I/Tf) versus required long term average uplink spectrum efficiency R m in/(NW), (b)-(d) I/Tf versus q for A = \i and various values of 

Rmin/(NW). 



collision time), caused by sensing error, is examined in the 
next subsection. 



C Simulation Results 

We present some simulation results in this subsection to 
demonstrate the effectiveness of the proposed ergodic CRN 
spectrum sharing strategy. The number of sub-channels is 16 
(TV = 16) and the number of ad-hoc bands is 4 (M = 4). Each 
ad-hoc band overlaps with 4 consecutive CRN sub-channels, 
and the 4 ad-hoc bands do not overlap with each other. The 
time fraction parameter in each frame a is set to 0.5. 

The channel coefficients h]^ (where i G and j G 

{r, d}, i 7^ j) are modeled as independent and identically 
distributed Rayleigh fading with zero mean and unit variance. 
We assume that the relay node is located right in the middle 
between the source node and the destination node. The large- 
scale path loss factor of all the wireless links is set to 4. 
Suppose that the signal-to-interference-plus-noise ratio (SINR) 



NWN* 



5dB. 



Cx- Th en, the 



of the source-destination link is set to 
Assuming that = N% and P^ ax 
SINRs of both the source-relay and relay-destination links are 

T . P*E{\h s > r \ 2 } P^K{\h^ d \ 2 } 1P7 , A . 

equal to ma fr^ 1 J = m ™ N wN* = 17dB according 
to the path-loss factor. The simulation results are obtained 
by averaging over 500 realizations of NSI (averaging over 
500 frames). According to our real-time spectrum sharing 
strategy in Section IV-B, the computation and signaling delay 
is negligible, which means we can choose 6 = 0. 

Figure 6(a) shows the simulation results of normalized 
long term average collision time (I/Tf) versus required long 
term average uplink spectrum efficiency R m i n /(NW), for 
XTf = fiTf — 1 for all m = 1, . . . , M. To compare with the 
proposed ergodic spectrum sharing strategy, we also perform 
the same simulations with its relay-free and sensing-free 
counterparts. The performance of the ergodic spectrum sharing 
strategy using only phase- 1 spectrum sensing is also presented. 
We can observe from Fig. 6(a) that the proposed strategy 



14 



10" 



10" 



10"' 











































































■rmT ___ n __ T _ n- 












Rmin/(NW) = 2.8 bits/s/Hz 

:: R m in/(NW) = 1.7 bits/s/Hz 

:: Rmin/(NW) = 0.6 bits/s/Hz 





































































10 



10 10 10" 

Sensing error probability 



10 u 



Fig. 7. Simulation results of normalized long term average collision 
time (I/Tf) versus sensing error probability, for different values of 

Rmin/(NW). 



outperforms both the relay-free strategy and the sensing- 
free strategy. Moreover, the proposed strategy with spectrum 
sensing in two phases performs better than with only phase- 1 
spectrum sensing. 

To further examine how the behavior of the ad-hoc traffic 
affects the performance of the proposed strategy, we define 
a parameter, called the relative variation rate of the ad-hoc 
traffic state or the relative sensing period of the CRN, as 



(53) 



where A = A and \i = \i for all m. A small value of q 
(that corresponds to small values of A and n) implies that 
the on-off state of the ad-hoc traffic changes slowly in each 
CRN frame. However, the ad-hoc traffic state would change 
many times in each CRN frame if q is large (that corresponds 
to large values of A and n). Figures 6(b) to 6(d) show the 
simulation results of normalized average collision time (I/Tf) 
versus relative variation rate q for A = /i and various values 
of R m in/(NW). Since the interference metric of sensing-free 
strategy is determined by that ratio A//i, but not how fast the 
ad-hoc traffic varies (see Eq. (37)), the normalized average 
collision time of the sensing-free strategy is constant versus 
q. From Fig. 6(b), one can observe that the proposed strategy 
performs best. However, the performance gaps between the 
proposed strategy and the relay-free and sensing-free strategies 
decrease with q, because the ad-hoc traffic is more difficult to 
predict for large q. For very large values of q, the proposed 
strategy has similar performance with the sensing-free strategy. 
Therefore, spectrum sensing provides no further benefit in this 
case. 

We can also see from Fig. 6(b) and Fig. 6(c) that the 
performance of the relay-free strategy seriously degrades as 
Rmin/(NW) increases from 0.6 bits/s/Hz to 1.7 bits/s/Hz. The 
performance degradation of the relay-free strategy is much 
larger than that of the proposed strategy because the CRN 
is capable of supporting higher uplink throughput than the 
relay-free strategy. In Fig. 6(d), the results of the relay- free 



strategy are not shown because this strategy is not feasible in 
supporting R m in/(NW) = 2.8 bits/s/Hz. 

We finally examine the robustness of our proposed strategy 
against sensing error. Under the same parameter setting associ- 
ated with the results shown in Fig. 6(a), Fig. 7 shows the sim- 
ulation results of normalized long term average collision time 
(I/Tf) versus sensing error probability, for different values of 
Rmin/ (NW). One can observe from Fig. 7 that, if the sensing 
error probability is small, e.g., less than 0.01, the performance 
degradation of the proposed strategy is insignificant. 

V. Concluding Remarks 

In this paper, we have investigated optimal spectrum sharing 
between cooperative relay and ad-hoc networks. Physical- 
layer resource allocation and MAC-layer spectrum access of 
the CRN are jointly optimized such that the average traffic 
collision time between the two networks is minimized while 
guaranteeing the CRN throughput requirement. Both frame- 
level design and ergodic design have been considered. For the 
latter design, a real-time implementation method of the optimal 
spectrum sharing strategy has been presented, by exploiting 
the structure of Lagrangian dual optimization solution. This 
implementation method has the following merits: 

1) Most computations are accomplished off-line, leaving 
only simple tasks for real-time computations. 

2) Although the sensing outcomes and channel gains are 
acquired at spatially separate nodes, the information 
exchange does not cause additional control delay. 

3) Almost all the computation loads at the source and 
relay nodes can be released, at minimal expense of an 
insignificant amount of information exchanges. 

4) Additional spectrum sensing in Phase 2 can be exploited 
to improve collision prediction. The relay node can queue 
up its received data if the sensing outcome in Phase 2 is 
unfavorable, which provides more flexibility for collision 
mitigation. 

Simulation results have been provided to examine the per- 
formance of the proposed strategy. We have found that good 
collision mitigation performance can be achieved if the ad-hoc 
traffic varies slowly and the required throughput of the relay 
network is not too high. The presented real-time implementa- 
tion techniques may also be useful for real-time transmission 
control of other wireless networks. 
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Appendix A 
Proof of Lemma 1 

For i G {1,2}, suppose that 9$ — 9$ = max{tf ,n G 
Mm} for some n' G Mm, i.e., sub-channel n' has the longest 
transmission time among the sub-channels in M m . Because 
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I^CU neA/m #,we have 

J Pr {X m (t) = l\X m (0) = x m } dt 

< I Pi{X m (t) = l\X m {0)=x m }dt. (A.1) 

If the following condition is satisfied 

4*) CI«, Vn e A/* m , (A.2) 

for m = 1, . . . , M and i G {1, 2}, then 1$ = U ney v m 
Then, equality holds in (A.l), and the interference is mini- 
mized. Therefore, J is determined by only 1$ . 
We further show that the condition 

IW =I#, VneA/* m , (A.3) 

is satisfied at the optimal solution to Problem (P), and 9$ = 
§m = max{9n\n £ Mm} for any n e M m - 

Suppose that (A.3) does not hold at the optimal solution to 
Problem (P). There must exist a sub-channel k G M m such 
that Of < 9$. Then, one can increase ef until ef = 9$ 
Ctt£° =0); this will increase Rcrn in (2) without changing 
the values of I in (6), because 1$ remains the same. In order 
to reduce /, one can further scale down §m to reduce / and 
Rcrn simultaneously, until equality holds for the constraint 
Rcrn > ^min- In summary, if (A.3) is not true, then one can 
always achieve a smaller objective value I for Problem (P). 
Thus, (A.3) is satisfied at the optimal solution to Problem (P), 
and the remaining problem is to determine the optimal 1^) . 

Following the arguments in [17, Lemma 1], one can further 
show that if x m = 0, the optimal transmission times are given 
by 4 1} = [ST f , (5 + 0$)T f ] and 1^ = [aT f , (a + 0$)T f ] 
for all n G M m , and if x m = 1, the optimal transmission 
time intervals are given by = [(a )!>,«!>] and 
ln ] = [(1 - 9$)T f ,T f ] for all n G M m - Lemma 1 is thus 
proved. ■ 
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